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SUMMARY 

Flow/damage surfaces can he defined in terms of stress, inelastic strain rate, and internal variables using a ther- 
modynamics framework. A macroscale definition relevant to thermodynamics and usable in an experimental pro- 
gram is employed to map out surfaces of constant inelastic power in various stress planes. The inelastic flow of a 
model silicon carbide/titanium composite system having rectangular, hexagonal, and square diagonal fiber packing 
arrays subjected to biaxial stresses is quantified by flow/damage surfaces that are determined numerically from 
micromechanics, using both finite element analysis and the generalized method of cells. Residual stresses from pro- 
cessing are explicitly included and damage in the form of fiber-matrix debonding under transverse tensile and/or 
shear loading is represented by a simple interface model. The influence of microstructural architecture is largest 
whenever fiber-matrix debonding is not an issue; for example in the presence of transverse compressive stresses. 
Additionally, as the fiber volume fraction increases, so does the effect of microstructural architecture. With regard to 
the micromechanics analysis, the overall inelastic flow predicted by the generalized method of cells is in excellent 
agreement with that predicted using a large number of displacement-based finite elements. 


I. INTRODUCTION 

The advent of man-made fiber-reinforced composite materials some forty years ago enabled the design of more 
efficient structures and greatly expanded the domain of engineered materials. Recent developments in the processing 
of unidirectional metal matrix composites (MMCs) provide new opportunities for engineers and materials scientists 
to tailor microstructural architecture for specific applications. For example, placing individual fibers in photo-etched 
grooves in foils of matrix material results in a very uniform microstructure. Since the grooves hold the fiber in place 
during consolidation, the microstructure can be engineered by simply specifying the foil thickness and groove pat- 
tern. Consequently, this paper addresses the effect that various microstructures have on the overall inelastic material 
response in the presence of multiaxial stress states. 

This research merges and builds upon two recent publications on unidirectional composites; one dealing with 
different microstructural architectures subjected to uniaxial loadings (Arnold et al., 1996a) and the other with 
macroscale flow/damage surfaces (Lissenden and Arnold, 1997a) given a fixed architecture. Hence, a brief summary 
of each follows, after which the objectives of the current work are enumerated. 

Arnold et al. ( 1 996a) begins with an extensive survey of the literature dealing with the effects of fiber shape and 
distribution on the response of composites. The generalized method of cells (Paley and Aboudi, 1992) was summa- 
rized and then used to predict the response of a silicon carbide/titanium system (SiC/Ti). The results presented dem- 
onstrate the relative effects of fiber arrangement and shape on the axial and transverse stress-strain and creep 
responses for both strong and weak fiber/matrix interfacial bonds. Fiber arrangements considered include rectangu- 
lar. square, square diagonal, and hexagonal periodic arrays as well as a random arrangement. Circular, square, and 
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cross-shaped fiber cross-sections were studied as were the effects of fiber volume fraction and strain rate on the 
stress-strain response. Additionally, the stress-strain and creep responses of “hybrid” composites (those with differ- 
ent size fibers having different arrangements and bond strengths) were investigated to determine the feasibility of 
using this approach to enhance the transverse toughness and creep resistance. 

Lissenden and Arnold (1997a) addressed theoretical and experimental issues regarding constitutive model 
development for anisotropic heterogeneous materials using micromechanics in lieu of actual experimental data. The 
generalized method of cells with four subcells, arranged to represent a square array of fibers was used to predict the 
multiaxial response of MMCs in terms of overall flow/damage surfaces (i.e., those incorporating matrix visco- 
plasticity and fiber/matrix debonding) in three different stress planes. The flow/damage surfaces were shown to have 
their centers offset from the origin by residual stresses and their shape altered by debonding. The normality condi- 
tion was shown to be reasonably well satisfied for initial surfaces of constant inelastic power (or in this case dissipa- 
tion) in the presence of fiber-matrix debonding. Results indicated which types of flow/damage surfaces should be 
characterized and what loading histories applied to obtain the most meaningful experimental data for guiding theo- 
retical model development and verification. This work differed from the extensive studies of plasticity in MMCs 
conducted by Dvorak and coworkers (summarized by Dvorak, 1991), in that different definitions of rate-dependent 
How and the effects of fiber-matrix debonding were addressed by Lissenden and Arnold ( 1997a). 

Herein we will demonstrate, using micromechanics, the influence that the continuously reinforced periodic 
microstructures shown in figure 1 (namely rectangular, hexagonal, and square diagonal fiber packing arrays) have 
on overall flow/damage surfaces in the axial-transverse (Gj | - o 02 ), transverse-transverse (G^ - O^), axial-shear 
(0| j - Op), and transverse-shear (o^ - Op) stress planes. Both a displacement-based finite element analysis (FEA) 
and the generalized method of cells (GMC) micromechanics approach are employed. Since the material response is 
in general viscoplastic, we will consider flow surfaces (typically important at elevated temperature), which are geo- 
metrically analogous to yield surfaces (typically important at room temperature). Furthermore, silicon carbide/ 
titanium (SiC/Ti) composites exhibit a weak bond between the constituents and this damage mode, when active, can 
greatly affect the overall material response. The effects of damage are included in the definition of flow surfaces, 
hence the name flow/damage surfaces. We will pay particular attention to definitions which indicate a prescribed 
level of inelasticity or damage, in an attempt to balance theoretical and experimental considerations. Sometimes 
results for the strong bond case (which is fictitious for current SiC/Ti sy stems) are shown in order to demonstrate 
how fiber/matrix debonding influences the predicted overall inelasticity. 

MMCs are anticipated to be typically used in elevated temperature environments. Unless stated otherwise, the 
results presented are for a TIMETAL-21S titanium matrix reinforced by a 35 percent fiber volume content of con- 
tinuous SCS-6 fibers. This unidirectional SiC/Ti system is taken as our model material system (see Arnold and 
Castelli, 1995). Two fairly extreme temperatures are considered; 23 and 650 °C (where time-dependent effects 
dominate the titanium matrix). The fiber response is taken to be linear elastic and temperature-independent, while 
the matrix response is elastic-viscoplastic (Arnold et al., 1996b and c) and highly temperature-dependent as illus- 
trated in figure 2. Note that the protective fiber coatings used to reduce fiber strength degradation (associated with 
chemical reactions with the matrix during processing) result in a weak bond to the matrix. A discussion of modeling 
this weak bond and the viscoplastic matrix response is postponed until Section III. Finally, as we are int erested in 
surfaces in prescribed stress planes, all loadings are stress controlled. An equivalent stress rate, ^Jc jd j > of 
2 MPa/sec is used in all cases, where Cj is the contracted stress tensor and / = 1,2,. .., 6 . 

In the next section (II) we review the theoretical framework and basic concepts applicable to inelastic flow' in 
MMCs. Flow surface definitions and factors that influence flow are discussed and an appendix (A) provides a con- 
tinuum example. The FEA and GMC micromechanics models are introduced along with the constituent constitutive 
models in Section III. The stress-strain response and inelastic power-time response for different microstructural ar- 
chitectures are presented in Section IV. A related appendix (B) discussts the accuracy of the FEA and GMC models. 
Finally, Section V presents the effects of microstructure and fiber-matrix debonding on macroscale flow/damage 
surfaces in four different stress planes at 23 and 650 °C. Section VI provides a closure for the paper. 


II. INELASTIC FLOW 

Our subject is the inelastic multiaxial response of MMCs for applications in elevated temperature environments. 
Hence, the viscoplasticity of anisotropic composite materials susceptible to internal damage must be addressed. 

From the viewpoint of a structural analyst investigating complex components, a macroscale continuum model is 
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preferred due to its numerical efficiency relative to alternative micromechanics analyses. However, such a con- 
tinuum model is not currently available; nor are the experimental results that are necessary to develop such a model. 
Thus, in the present study we will utilize a theoretical framework appropriate for continuum modeling to guide us 
and employ two micromechanics models to simulate the required experiments numerically, so as to understand the 
theoretical and experimental implications and/or assumptions necessary in either a macro (continuum) or 
micromechanics based approach. 


II. 1. Theoretical Framework 

The theoretical considerations used in this paper are based on an energy balance. The primary variables at any 
point are the Cauchy stress tensor, G-, the internal stress tensor, a y/ , and temperature, T. Other internal state vari- 
ables could also be defined and used if we were so motivated. The current values of these variables can be used to 
define the Gibbs thermodynamic potential, G - Gfo^a^D^T), where D f/ is the preferred direction tensor formed by 
the self product of the unit vector denoting the fiber direction. Conjugate to these variables are the total strain tensor, 
t ij% the internal strain tensor, and the entropy, S, 


dG dG 

£ij ~ dOjj ' ij ~ doty' ~ dr 


( 1 ) 


Our basis is that the total work performed on the system must be equal to the sum of the stored energy and the en- 
ergy dissipated, where the stored energy includes an elastic component as well as an inelastic component associated 
with the internal state. Thus, the dissipation potential, Q = QiG^a^D^T), can be defined to be GyE-ij — CLyAjj. 

The associated How law (for the inelastic strain rate) and evolution equations (for the internal strain rate) are 
given by normality. 


and 


£ h 


aa 

do,/ 


; _ dn 

* ij da , 


tty - Qijkl^kh 


( 2 ) 


where 


Q ‘i U ‘ daqtou 

and Qj k( is called the internal compliance operator (Arnold and Saleeb, 1994). Thus, once the functional dependen- 
cies o/ the Gibbs and dissipation potentials have been determined, all of the variables are known by simple 
differentiation. Let us now assume that the dissipation potential can be written in terms of two scalar functions, 

Q = £2(/\/f), where F depends on the deviatoric effective stress, Z-, and H depends only on the internal stress 
(Robinson and Ellis, 1986). The deviatoric effective stress is the difference between the deviatoric Cauchy stress and 
the deviatoric internal stress. Now the How law can be written. 


. / _ 3Q dF 
£,> “ 3 F do,, ' 


(3) 


Thus, the direction of the inelastic strain rate vector is normal to surfaces having F = constant. However, if it is 
true that the normality condition is not satisfied in MMCs, as indicated by Nigam et al. (1994a,b) for boron/alumi- 
num, then it becomes necessary to develop a nonassociated How law and evolution equations. We return to this issue 
at the end of Appendix A. 
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II.2. Surface Definitions 


The concept of a yield surface is well known in rate-independent p asticity, even if no one definition of yielding 
has been universally adopted. The most common definitions employed are the proportional limit, a small (usually 
5 to 20 pm/m) offset strain, a back-extrapolation, and a large (usually 2000 |im/m) offset strain. For rate-dependent 
plasticity (viscoplasticity), the concept of a strict yield surface breaks down as stress states outside the yield surface 
are accessible (since no consistency condition applies). Thus the need for geometrically analogous, thermodynami- 
cally based, flow surface definitions. At least two different rate-dependent definitions have been proposed for mono- 
lithic materials (Clinard and Lacombe, 1988); 


( 1 ) surfaces of constant dissipation rate (SCDRs), defined by OjjEj/ — O-ij^ij and 

(2) surfaces of constant inelastic strain rate (SCISRs), defined by sjjtjj . 

These definitions can also be applied to composite materials. However, one needs to be careful when using 
micromechanics as the inelastic strain is an eigenstrain and consequently; 


<J, * — f OijijjdV 

l J l J y J y l J l J 

(4) 

ijjifi * — f e/.e/. JV 

v v y J y l J l J 


where an over-bar indicates a macroscale quantity and V is the total volume of the representative volume element. 
The difference is due to the existence of residual stresses, at least some of which are associated with stored elastic 
energy. Being an additive quantity, the overall dissipation potential can be taken to be the volume average of the 
local dissipation potentials (Suquet, 1987), and from a theoretical standpoint that is what we would like to quantify. 
However, it is indeterminate experimentally because the local (internal) variables are unknown and not measurable. 
This issue is dealt with more completely by Lissenden et al. (1998). Here let us define another type of flow surface 
using experimentally measurable macroscopic quantities, o f y£^, that are available from both micromechanics and 
macroscale models. We will call them surfaces of constant inelastic power (SCIPs). This is in fact the nomenclature 
that should have been used by Lissenden and Arnold (1997a) instead of the term initial SCDRs. The actual differ- 
ence between SCIPs and SCDRs depends on many factors, but as shown in figure 3 for transverse tensile loading, it 
can be substantial. 

Lissenden and Arnold (1997a) demonstrated, using micromechanics, that the direction of the overall inelastic 
strain rate vector can differ significantly from the outward normal of a SCISR. Whereas, for the stress planes consid- 
ered, the direction of the overall inelastic strain rate vector was reasonably close to the outward normal of the 
SCDRs (actually SCIPs) 1 considered. SCISRs, however, are certainly more amenable to experimental methods than 
are SCDRs, and may be more amenable than SCIPs given that stress quantities are not included in the definition. For 
macroscale continuum theories it is not possible to distinguish between SCDRs and SCIPs. A macroscale continuum 
example is given in Appendix A. 


II. 3. Factors Influencing Flew 

Many factors influence inelastic flow in metallic materials. Certain y, temperature and loading rate as well as 
the past loading history can be important in many metals. Additionally, microstructural architecture, degree of 


'Note that for a near virgin state, SCDRs and SCIPs are nearly the same given a- ~ 0, bu? as incleastic flow occurs (tt - * 0) they will diverge. 
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anisotropy, fiber-matrix bond strength, and damage influence flow in composites. Other factors, such as the stress 
plane, definition, and target value influence how flow is represented. Lissenden and Arnold ( 1997b) illustrated the 
effects of many of these factors. The current paper focuses on the influence of microstructural architecture and fiber- 
matrix bond strength. We consider repeating microstructures, specifically ones having rectangular, hexagonal, and 
square diagonal fiber packing arrays as shown in figure 1 . The rectangular array has an aspect ratio, R = alb. For the 
special case of a square array, R = 1. Additionally, we consider strong and weak liber-matrix bonding. Our 
definition of strong is that there is no discontinuity in the displacement field at the interlace. Likewise, weak means 
that the interface can transmit a finite traction before debonding causes the fiber and matrix to separate. 


111. MICROMECHANICS 

Biaxial experiments on unidirectional continuous-liber reinforced SiC/Ti in the axial-transverse (0| | — CJ 22 ), 
transverse-transverse (o-,-,- o 33 ), axial-shear (o, , - 0 12 ), and transverse-shear (0 22 -0| 2 ) stress planes were simu- 
lated numerically using micromechanics. Initial overall (maeroscale) flow/damage surfaces were mapped out by a 
sequence of stress-controlled proportional loading probes at different angles in a prescribed stress plane as shown in 
figure 4. The two nonzero stress rate components were calculated to be the cosine and sine ol the equivalent stress 
rate (2 MPa/sec). Each probe started at the origin and continued until the inelastic power target value had been 
reached, after which the material was returned to its virgin state, the probe angle increased (usually by 5°), and the 
next probe conducted; until the entire stress plane (0 -> 360°) had been probed. For surface determinations at room 
temperature thermal residual stresses were accounted lor by cooling the composite from a stress-lree temperature of 
815 °C in 2 hr. These simulations were conducted using both FEA and GMC, with the specifics described below. 


III. I . Finite Element Analysis 

The commercial FEA program ABAQUS (HKS.1995) was used to determine the response of the microstruc- 
tures of interest (fig. 1 ). Each microstructure can be represented by a repeating unit cell that has two planes of sym- 
metry. Thus, it is only necessary to analyze a quarter of each unit cell. Coarse and fine discretizations of one quarter 
of the repeating unit cell for both square and square diagonal arrays are shown in figure 5. Only strongly bonded 
SiC/Ti was considered using FEA for determining SCIPs. Generalized plane strain triangular elements formulated 
by Lissenden and Herakovich (1995), and implemented into ABAQUS through a UEL subroutine, were used for 
modeling the repeating unit cell. These triangular elements have an extra node that is common to all elements in the 
mesh. The degree of freedom in the 1 -direction (axial) at this common node is what makes this a generalized plane 
strain element. Furthermore, warping of the 23-plane permits axial-shear loading to be simulated. Overall inelastic 
strain components were calculated as the difference between the total strain components, found by volumetric aver- 
aging, and the elastic strain, found from the overall stress (volumetric average) and elastic properties. 

The boundary conditions applied to each FEA mesh are illustrated schematically in figure 6. To maintain com- 
patibility with the adjacent unit cells, multipoint constraints (MPCs) were applied on the edges that are in contact 
with adjacent unit cells. For example, MPCs were applied along the right edge such that the edge remains straight. 
Stress control loading was simulated by applying concentrated forces at the nodes along an appropriate edge; except 
for axial loading, where a concentrated force was applied to the common node, which is shown in figure 6. Details 
of the applied loadings are summarized in table I. 


II1.2. Generalized Method of Cells 

The generalized method of cells (GMC) (Paley and Aboudi. 1992; Aboudi, 1995) is an approximate analytical 
micromechanics model that extends the original method of cells (Aboudi, 1991) to an arbitrary number ot subcells, 
permitting the study of different microstructures. The reader is referred to Aboudi (1995) for the derivation of the 
detailed equations of GMC. The primary equations of GMC relate to standard micromechanics equations as follows 
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e ij (( e //0)))’ 

ot=l P=l y= 1 

(5) 

u=l (5=1 7 =| 

(6) 

CT i) = “£»)• a ij~c 

(7) 

(CW*VW*>)). c * =^X X J i d a h p l r c'«l , y>A^r< 

a= 1 P=\ /=! 

(8) 



= -S (C * r 'X X 

a=l p= 1 7=1 

(9) 

*!j = {{%;/(*)£«(■.))). 


(X=\ p=\ 7=1 

(10) 


In equations (5) to ( 10) an overbar indicates a macroscale variable; (x) explicitly indicates that the local variable is a 
function of position; {{■)) = — J ■ dVis the volume average operator; C * A7 denotes the overall effective elastic stiff- 
ness tensor; the superscript (a|Jy) refers to subcell (apy) where a, p, and y are indices in the three coordinate direc- 
tions; d h , and / are the overall dimensions of the repeating unit cell and the same variables with subscripts denote 
dimensions of subcells; A ijkj> B ijkl , and D y kl are the elastic strain, stress, and nonelastic strain localization tensors; 
and e/and e/ are arrays containing the subcell inelastic and thermal str; in vectors. 

GMC has been implemented into the recently developed micromechanics analysis code (MAC/GMC) which 
has many user friendly features and significant flexibility (Wilt and Arnold, 1996). The GMC unit cell representa- 
tions utilized herein for square, hexagonal, and square diagonal arrays are shown in figure 7. 


III. 3. Constituent Models 

As mentioned previously, the fiber response is assumed to be linear elastic and temperature independent 
(E — 400 GPa, v = 0.2, CTE = 3.5x1 °C -1 ). The elastic-viscoplastic tehavior of the matrix is represented using a 
generalized viscoplastic potential structure (GVIPS) model (Arnold et a ., 1996b,c). This model is a fully associa- 
tive, multiaxial, nonisothermal, nonlinear kinematic hardening viscoplastic model for use with initially isotropic 
metallic materials. A unique aspect of this model is the inclusion of nonlinear hardening through the use of a com- 
pliance operator Q jjk! in the evolution law lor the back stress. This nonlinear tensorial operator is significant in that it 
allows both the flow and evolutionary laws to be fully associative and greatly influences the multiaxial response 
under nonproportional loading paths so prevalent w ithin composite mate rials. The primary equations are summa- 
rized by. 
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flow law 


where 


e'j = XL i} if F > 0 

_ [Qijkl^-klpqb pq *f 

'M b,j ^ 

a tj = Lijklftu -Okl?) 


a ijZij <0 


evolution law 

internal constitutive rate law 



(16) 


where the Macau ley brackets are defined by 


The constants K , - K 4 contain material parameters and internal stress invariants and can be interred from Arnold 
et ah (1996b). 

Weak bonding between the fiber and matrix is modeled by assuming that a jump in the displacement field may 
occur under certain conditions, while the traction vector remains continuous. In this model debonding initiates when 
the normal traction exceeds a critical value or when the tangential traction exceeds a critical value, with no interac- 
tion between the two. 


«„ = R ,j„\ 
>*, = R j,\ 


if 


T > T n "or 

l n — * t\ 


T t > Tj C 


(17) 


where u denotes interfacial displacement, T the interfacial traction, and R the flexibility of the failed interface. The 
subscripts n and t denote the directions normal and tangent to the interface respectively. This simple model com- 
bines elements of interface models developed by Jones and Whittier (1967) and Achenbach and Zhu (1990). Once 
debonding has initiated, the interfacial displacement rate is made proportional to the stress rate, simulating a per- 
fectly plastic-like behavior. In the results for SiC/Ti with a weak bond the critical normal and tangential tractions 
have been taken to be 103 and 41 .4 MPa, respectively, and the flexibility of the failed interface has been taken to be 
0.271 mm/MPa. 


IV. STRESS-STRAIN RESPONSE 

Before determining flow surfaces we validated our FEA models by performing a convergence study and com- 
paring stress-strain responses with those published by Brockenbrough et al. (1991) for boron/aluminum. In figure 8, 
the stress-strain and inelastic power-time responses are shown for axial (Oj |), transverse (C^)* an< ^ shear load- 
ings at 650 °C. Here the effect of microstructure, i.e., square, hexagonal, and square diagonal repeating arrays of 
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fibers, is demonstrated using CMC and the unit cells shown in Figure 7. Both strong and weak fiber-matrix bonds are 
considered, but as expected for axial loading the interface does not affect the response. Also as expected, micro- 
structure has no affect on the axial response (fig. 8(a)). Hexagonal and square diagonal array responses for trans- 
verse and shear loadings (figs. 8(b) and (c)) for strongly bonded composites are very similar and diverge from the 
response of a square array for transverse loading, but agree closely with the response of a square array for shear 
loading. It is clear that the transverse response is strongly affected by microstructural architecture, while the axial 
and shear responses are not affected and very mildly affected, respectively. However, for weakly bonded composites 
the square and hexagonal array transverse and shear responses are similar, while the square diagonal array exhibits a 
more compliant response (figs. 8(b) and (c)). 

An accuracy comparison between these CMC unit cell representations and their FEA counter-parts was also 
undertaken and is described in Appendix B. A summary of our observations based on stress-strain and inelastic 
power-time responses, relative to the appropriate micromechanics model (FEA and GMC) and microstructural archi- 
tecture (square, hexagonal, and square diagonal arrays) is as follows: 

• GMC predictions agree well with FEA predictions if a fine FEA mesh is employed; 

• A fine FEA mesh is required to accurately predict gross inelastic flow in the transverse direction, but for axial 
and shear loadings a coarse FEA mesh is adequate, however, a coarse mesh appears adequate for any of the 
loadings considered, provided the inelastic power is small ; 

• Microstructural architecture is very significant for transverse loading, but has much less effect for axial 
(none) and shear loading (only significant for weakly bonded composites). Note that we only considered axial 
shear (Op) loading, not transverse shear (G^) loading. 


V. MACROSCALE FLOW/DAMAGE SURFACES 

Our focus now turns to macroscale tlow/damage surfaces in order to obtain information on the initial overall 
dissipation potential. Surfaces of constant inelastic power (SCIPs), defined simply by will be determined in 

order to give this work relevance to experimental and macroscale continuum approaches. Here the inelastic strain 
tensor includes matrix inelasticity as well as nonlinear effects associated with fiber-matrix debonding. First, we 
exercise the FEA models for a strongly bonded composite, then compare FEA and GMC predicted flow surfaces, 
and then finally consider a weakly bonded composite using GMC. 

Consider first the effect that the target value has on SCIPs in the axial-transverse stress plane at 23 °C. Figure 9 
shows 1, 5, and 10 kPa/sec SCIPs for both a square array and a square diagonal array as predicted by the fine and 
coarse FEA meshes shown in figure 5 for strong bonding. The first observation that we make is that the surfaces for 
different target values are not concentric for either array and that they i re offset in the axial compression direction 
due to thermal residual stresses. While the surfaces for square and square diagonal arrays are approximately the 
same size for the same target value, their shapes are different, particularly the 5 kPa/sec SCIPs in the first quadrant 
(i.e.. tension-tension). The 5 kPa/sec SCIP for a square array has a 'nose’ at approximately 35°, while the ‘nose' on 
the same SCIP for a square diagonal array is at approximately 15°. The differences in shapes of the overall flow' 
surfaces are apparently related to the different local stress and strain fields. 

Local J~> (1/2 SySp where S {J is the deviatoric Cauchy stress) contours in the matrix are plotted in figure 10 for 
both arrays given their fine mesh idealizations and all three target valies for loading at 35° from the axial stress 
axis. The smallest values (darkest) occur in the matrix adjacent to the iber near the bottom edge (and top edge for 
the square diagonal array). This is interesting because 7 ? is small whete the transverse stress is large and large 
where the transverse stress is small. Another way to see the differences in the local fields for square and square 
diagonal arrays is to compare the local stress invariants / j (the sum of the normal stress components) and 7-,. In 
figure 1 1, the stress invariants / j and J-, for each matrix element have been sorted and are plotted in descending 
order. Note that while the two arrays each have a fiber volume fraction of 0.35, the square array has 227 matrix 
elements in the fine mesh, while the square diagonal array has 281 . The square diagonal array has a relatively 
small number of elements with large hydrostatic stress, while the square array has a more uniform distribution of 
hydrostatic stress (fig. 1 1(a)). For example, the difference between tht square and square diagonal arrays for 
5 kPa/sec SCIPs is related to the different /, distributions shown in figure 1 1(b). The square array has fewer ele- 
ments with high values of/, and therefore requires a higher overall stress to achieve the same How target value as 
the square diagonal array. 
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Comparing the SCIPs from the fine and coarse mesh FEAs shown in figure 9 indicates that the coarse mesh for 
a square array of fibers provides reasonably good results in the axial-transverse stress plane. However, in the case of 
a square diagonal array of fibers agreement between the fine and coarse mesh results is not as good, implying that 
too few elements were used in the coarse square diagonal array mesh. 

Clearly, it is highly preferable to use a coarse mesh in lieu of a fine mesh, provided the accuracy of the results 
can be maintained, as the execution time is greatly reduced. 

Nimmer et al. ( 1991 ) used FEA to study the effect that the aspect ratio of a rectangular array has on the trans- 
verse tensile response of a weakly bonded SCS-6/Ti-6-4 composite system. Here we expand the spirit of Nimmer s 
study to the biaxial loading behavior of a strongly bonded SCS-6/TI METAL 21 S composite system. A wide range 
of aspect ratios (R - a/b in fig. 1 ), 0.5 < R < 2.0, are considered and compared with results for a square array 
(R = 1.0). Figure 12 shows 5 kPa/sec SCIPs at 23 °C in four stress planes for aspect ratios of 0.5, 1 .0, and 2.0. As 
demonstrated previously, the microstructural architecture, defined by the aspect ratio here, has the largest effect by 
far on the transverse response. 

One application for the kind of information that How surfaces provide is the design of components subjected to 
loads resulting in deterministic multiaxial stress states. The microstructure can then be engineered to best resist the 
loading. For example, a ring mounted on a shaft in a jet engine will rotate in service. Thus, the radial and circumfer- 
ential stress components will be tensile. Suppose the ring is to be fabricated from hoop- wound, strongly bonded SiC/ 
Ti, and that the goal is to delay the onset of inelastic How as long as possible. According to the first quadrant of the 
axial-transverse stress plane in figure 12, the largest possible aspect ratio should be used. However, for a weakly 
bonded composite, once the interface fails the stress distribution is much different because the transverse stress in 
the fiber gets redistributed to the matrix. In this case, a small aspect ratio is preferred because it provides a long liga- 
ment of matrix material between fibers for the transverse stress to flow through, while a large aspect ratio is associ- 
ated with a short ligament of matrix between fibers and more localized flow (Lissenden and Herakovich, 1996). 

In figure 13, FEA (coarse mesh) and GMC predicted 1 and 10 kPa/sec SCIPs at 23 °C are compared in all four stress 
planes for composites with a strong fiber-matrix bond. The SCIPs predicted by FEA and GMC are in good agree- 
ment. This is in line with the excellent agreement between yield surfaces defined by local yielding (Mises stress) 
that Pindera and Aboudi (1988) reported for the method of cells (square packing) and FEA. In figure 13 overall 
SCIPs are determined based on overall inelastic strain rates and stresses, as might be done in an experiment. These 
results are in contrast to the FEA/GMC comparisons made in Lissenden and Arnold (1997c), where GMC and FEA- 
predicted flow surfaces were obtained at 650 °C using strain-controlled and stress-controlled loading, respectively. 
While initial results indicated that How surfaces in the axial-transverse stress plane were not very sensitive to 
whether stress or strain control probing was used (Lissenden et al., 1998); in the transverse-transverse stress plane 
the control mode is important. Thus, the larger flow surfaces obtained by GMC in strain-control (fig. 8 in Lissenden 
and Arnold, 1997c) are easily explained. Inelastic deformation is strain rate dependent, thus more deformation oc- 
curs for a slower overall strain rate than a faster one. Consequently, in strain control the overall strain rate is fixed, it 
is the stress rate that decreases as inelastic deformation occurs. However, in stress control the overall strain rate must 
decrease as inelastic deformation occurs, resulting in more inelastic flow relative to strain control. 

Figure 14 shows excellent agreement between GMC and FEA (coarse mesh) 10 kPa/sec SCIPs at 650 °C in all 
four stress planes. The sharp points apparent on some of the surfaces are artifacts due to the large step in probe angle 
(5°) used to make the analysis less time consuming. These SCIPs are actually smooth curves. SCIPs from finer 
GMC discretizations having 16 subcells (a cross shaped fiber) and 49 subcells (a roughly circular fiber) were also 
examined in the axial-transverse stress plane. As was found for FEA (fig. B 1 ), the axial response was independent of 
the discretization, but the transverse stress was 7 percent less for the 49 subcell model than for the 4 subcell model 
shown in figure 14. This difference is almost imperceptible at the scale it is drawn. Local stress and strain fields, 
which ultimately determine the strength or cyclic life of a material, are not reported here, but are the topic of a com- 
panion paper (Lissenden et al., 1998). 

Consider now, the case of a weak fiber-matrix bond using GMC, which is more realistic for the current SiC/Ti 
system. In figure 15, 1 kPa/sec SCIPs at 23 °C for a square array are shown for weak and strong bonds. The primary 
effect of a weak bond is to significantly reduce the tensile stress at which deviation from proportional response 
begins. Since compressive interfacial tractions are not detrimental to the integrity of the interface, debonding does 
not occur for transverse compressive loading (unless it is due to Poisson expansion, see Lissenden and Arnold 
(1997a)). Additionally, axial loading is not observed to cause debonding. 

The effect that microstructural architecture has on 1 kPa/sec SCIPs in the axial-transverse and transverse-shear 
stress planes for weakly bonded SiC/Ti al 650 °C is shown in figure 16(a). Microstructural architecture has the most 
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influence when a compressive transverse stress is present. As was obseved for the stress-strain response (fig. 8), the 
square array exhibits less inelastic flow than the hexagonal and square diagonal arrays. 

Brockenbrough et al. ( 1 99 1 ) and Arnold et al. ( 1 996a) demonstrated that as the fiber volume fraction increases, 
the effect of microstructural architecture increases for uniaxial loadings. Figure 16(b) shows 1 kPa/sec SCIPs in the 
axial-transverse stress plane at 650 °C for square, hexagonal, and square diagonal packings and a fiber volume frac- 
tion of 0.50. The primary effect of increasing the fiber volume fraction is to enlarge the SCIP. Additionally, larger 
differences in the SCIPs for the three microstructures are observed for the higher fiber volume fraction, especially 
where the transverse stress is compressive (and the interface has not debonded). 


VI. CLOSURE 

Micromechanics has been used to generate overall flow/damage surfaces for unidirectional MMCs having dif- 
ferent microstructures. The generalized method of cells results compared very well with those from the comput- 
ationally more demanding finite element analysis approach, wherein many elements are required to obtain the same 
level of accuracy. Except for the special case of uniaxial loading in the direction of the fibers, microstructure influ- 
ences the shape and location of the macroscale flow/damage surface. The magnitude of the effect that microstructure 
has on flow/damage surfaces depends strongly on fiber/matrix bond strength, fiber volume fraction, and the type of 
loading, among other secondary factors. A weak fiber/matrix bond significantly influences the shape of the flow/ 
damage surface. A dramatic flattening of the surface is observed for transverse tensile loading. Debonding also 
occurs under shear loadings, but its effect is less drastic. 

Flow/damage surfaces are important because multiaxial stress states are encountered in the analysis and design 
of most structural components due to complex loadings and geometries, as well as due to stress concentrations. The 
response of MMCs to multiaxial stress states can not be inferred from uniaxial test data. It is necessary to conduct 
multiaxial tests to guide and later validate theoretical models. To implement the thermodynamically based frame- 
work discussed in Section II for a macroscale continuum model, the functional form for the Gibbs and dissipation 
potentials must be known. Neither of these potentials can be directly determined experimentally. However, one pos- 
sible approach for continuum model development is to experimentally determine SCIPs in various stress planes and 
use a micromechanics model to interpret the results. The functional form of the Gibbs and dissipation potentials 
could then be inferred from the micromechanics model having an inelastic power that correlates well with the 
experimental results over a wide range of conditions. This work has shown that in addition to the dependence of 
these potential functions on stress invariants, internal state invariants, and temperature, as for isotropic monolithic 
materials, they may also depend on directionality, microstructure, and bond strength for MMCs. This adds yet an- 
other level of complexity to the macro-based continuum approach for modeling multiaxial inelastic response of 
MMCs. 
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APPENDIX A 


A MACROSCALE CONTINUUM EXAMPLE 

Consider, a unidirectional MMC and smear the microstructure such that we have a homogeneous anisotropic 
material. Thus, all variables in this appendix refer to composite quantities, there are no local (micro) variables. The 
effect of hydrostatic stress on the flow of most metals is quite small and usually neglected, as is the case here, even 
though composites can be dependent upon hydrostatic stresses (e.g., Aboudi, 1991 and Jansson, 1995). Assume that 
the dissipation potential can be written in terms of two scalar functions. Q. = Q( /%//), where H is a function of the 
deviatoric internal stress, and the reinforcement direction, defined by /U. and 
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and k T is the transverse shear strength; also. T| and CD denote threshold stress ratios for shear and normal loading, 
respectively. The constant c scales the effect of ./, on inelastic flow and must be determined experimentally. This 
form of F is for transversely isotropic materials and is due to Robinson and Dully (1990), except here Jj ^ as 
replaced with Drucker’s ( 1949) J~>/J 3 form (c = -9/4). In this theory, the inelastic strain rate vector is normal to sur- 
faces having constant F in stress space (see eq. (3)). Thus, normality only applies to flow surfaces proportional to 
surfaces having F = constant. It can be shown that 
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in the neighborhood of the virgin state, which are proportional to F = constant surlaces. On the other hand, SCISRs 
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For an isotropic material (that is, r| = to = 1 ), SCISRs reduce to 


TT _ V 2 j-, dn 

£,]£,i ~ k\F~+\)dF 


(F+\) + 


cJi 


/t 6 (F+!)H 27 


l + ^c 


1/2 


(A4) 


and if the J 3 contribution to inelastic flow is negligible (that is, c - 0) we get 


TT _ a/2(F+1) <Kl 


£■•£■ — 
IJtJ 


dF 


(A5) 


Thus, SCISRs are only proportional to surfaces having F - constant for the special case of isotropic A materials, 
and therefore the normality criterion only applies to SCISRs for this special class of materials. Also, Lissenden and 

/ [ 2 7 7 

Arnold (1997a) extrapolated these results to rate-independent plasticity by considering <T yfy and J— , where 

the latter is the usual equivalent inelastic strain definition often used to define a yield surface, and found results 
analogous to those for SCDRs and SCISRs. Hence, the Nigam et al. (I994a,b) results on boron/aluminum that the 
plastic strain rate vector was not normal to the yield surface could be merely dependent on their definition of yield. 
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APPENDIX B 


A MESH REFINEMENT STUDY 

Here a comparison of the stress-strain responses for uniaxial and biaxial loadings using both FEA and GMC is 
conducted. Additionally, since we are interested in inelastic power, its increase as a function of time will also be 
plotted for these loadings. In figure B 1 the stress-strain and inelastic power-time responses are shown for axial 
(a, | ), transverse (G ?0 ), and shear (o p ) loadings at 650 °C. A square array of fibers is discretized using the coarse 
and fine FEA meshes (fig. 5) and the four subcell GMC model (fig. 7). For axial loading the three predictions are 
nearly identical for both stress-strain and inelastic power (fig. Bl(a)), which is not surprising since the material 
response is dominated by the fiber as evidenced by the high stresses and small inelastic power relative to the other 
two types of loading. On the other hand, transverse and shear loadings are matrix dominated and differences be- 
tween analysis techniques are anticipated. 

The disparity in the coarse and fine mesh FEA-predicted transverse stress-strain and inelastic power-time 
responses shown in figure Bl(b) indicates that the coarse mesh is not adequate for predicting gross inelastic flow, 
that is, the coarse mesh has not yet converged. Further refinement of the fine mesh indicated that it had converged 
for the range of stresses obtained in the surface determinations presented in this paper, but for larger stresses an even 
more refined mesh would be necessary. The GMC-predicted responses agree very well with the fine mesh FEA pre- 
dictions for the range of stresses obtained in surface determinations. However, for larger stresses a more accurate 
geometric representation of the fiber within the repeating unit cell would provide more accurate results, see figure 1 
of Arnold et al. 1996a. 

The shear responses (fig. B 1 (c)) from the two FEA meshes are identical, indicating that the coarse mesh has 
converged for shear loading. The GMC predicted response exhibits more inelastic flow and is in reasonably good 
agreement with the FEA results. The exponential form of the inelastic power accumulation with time (fig. B 1 (c)) is 
a result of the constant stress rate used, which causes the inelastic strain rate to approach infinity as the slope of the 
stress-strain curve approaches zero. 

Stress-strain and inelastic power-time responses calculated utilizing FEA and GMC are shown in figure B2 for 
proportional biaxial loading of a square array of fibers at 650 °C. The axial-transverse (G j , = G 22 ) response shown 
in figure B2(a) indicates that each micromechanics approach is equally good for predicting axial response, but that 
differences are present in the predicted transverse response, which lead to differences in the inelastic power. As for 
uniaxial transverse loading, the coarse FEA mesh is inadequate for predicting gross inelastic How while GMC pro- 
vides a very reasonable approximation of the inelastic flow. Each model is equally good for predicting the trans- 
verse-transverse (o-n = G 33 ) response, as shown in figure B 2 (b) and reasonable agreement is also observed lor the 
axial-shear (Gj , = g 12 ) response shown in figure B2(c). 

Note that the finely meshed FEA takes a relatively long time to execute compared to that of the coarse FEA 
mesh (a clock-time ratio of approximately 4: 1 ), and even more so relative to GMC. This lengthy execution time is 
greatly exacerbated when performing flow surface determinations as at least 72 directions for one surface determina- 
tion are required. Fortunately, the coarse FEA mesh is reasonably accurate for small values of inelastic power, 
approximately 10 kPa/sec and below. Thus, it appears acceptable to use the coarse mesh for surface determinations 
using transverse stresses provided the target value is small. Additionally, GMC, which is numerically more efficient 
than either FEA model, agrees well with the fine meshed FEA results. 
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Figure 1 . — Periodic microstrctures; (a) rectangular, (b) hexagonal, and 
(c) square diagonal array. 



Figure 2. — Predicted constituent response for SiC 
fiber and Ti matrix. 
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SiC/Ti at 650°C 
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Figure 3. — Inelastic power and dissipation for transverse 
tensile loading. 
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Square array (coarse mesh) Square diagonal array (coarse mesh) 




Square diagonal array (fine mesh) 


Figure 5. — Coarse and fine mesh FEA discretizations of the rectanglar, hexagonal, 
and square diagonal microstructures. 
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Figure 6. — Boundary conditions for FEA. 



Figure 7. — GMC discretizations of the rectanglar, hexagonal, arid square diagonal microstructures. 
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Figure 8. — Uniaxial stress-strain and inelastic power-time responses at 650 °C for square, hexagonal, and 
square diagonal arrays from GMC. 
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Transverse Stress (MPa) 



Axial Stress (MPa) Axial Stress (MPa) 


Figure 9.— Effect of target value on SCIPs at 23 °C in the axial -transverse stress plane for square and 
square diagonal arrays discretized with fine (dark lines and o) and coarse (light lines and □) FEA meshes. 
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Figure 14. — Comparison of 10 kPa/sec SCIPs pre 
for a square array. 
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Figure 16. — Effect of microstructure and fiber volume fractron on 1 kPa/sec SCIPs 
predicted by GMC at 650 °C given weak fiber-marix interface. 
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Figure B1 . — Uniaxial stress-strain and inelastic power-time responses at 650 °C for a 
square array using a coarse mesh FEA, a fine mesh FEA, and GMC. 
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Figure B2. — Biaxial stress-strain and inelastic power-time responses at 650 °C for a 
square array using a coarse mesh FEA, a fine mesh FEA, and GMC. 
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